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Abstract 

The NLO cross section of the gluon fusion process is matched to parton showers in the 
MC@NLO approach. We work in the framework of MadGraph5_aMC@NL0 and document the 
inclusion of the full quark-mass dependence in the SM as well as the state-of-the-art squark 
and gluino effects within the MSSM embodied in the program SusHi. The combination of the 
two programs is realized by a script which is publicly available and whose usage is detailed. 
We discuss the input cards and the relevant parameter switches. One of our focuses is on the 
shower scale which is specifically important for gluon-induced Higgs production, particularly 
in models with enhanced Higgs-bottom Yukawa coupling. 



1 Introduction 


Higgs production proceeds predominantly through gluon fusion in a large number of theories, 
including the Standard Model (SM). The recently discovered resonance in searches for a 
Higgs boson is fully consistent with the SM pictur^ so far. Still, the measured Higgs boson 
may be embedded in an enlarged Higgs sector with respect to the one of the SM which predicts 
only a single physical particle breaking the electro-weak symmetry. Two-Higgs Doublet Models 
(2HDM’s) such as the Minimal Supersymmetric SM (MSSM) are among the most popular theories 
with enlarged Higgs sectors. Such theories inevitably require the existence of further physical 
Higgs particles. A 2HDM predicts three neutral Higgs bosons: a light (h) and a heavy (H) scalar, 
and a psendo-scalar (A); as well as two charged Higgs particles (H^). Almost all 2HDM and 
MSSM scenarios that are in agreement with the experimental bounds feature a light scalar which 
is SM-like in its couplings to vector bosons and fermions, while the other Higgs bosons are heavier 
and, therefore, escaped detection up to now. Indeed, the experimental search for other Higgs 
resonances is one of the major focuses regarding the discovery of physics beyond the SM (BSM) 
in the second run of the Large Hadron Collider (LHC). 


Higgs production through gluon fusion is mediated by a colored particle. In the SM, the top quark 
gives the dominant contribution to the cross section [6]-[^. While also the bottom quark gives a 
sizable contribution, the effects due to other quarks are small and therefore usually neglected. In 
the 2HDM and the MSSM the Higgs-bottom Yukawa coupling can be enhanced with respect to the 
one of the top quark and the bottom loop may even constitute the dominant contribution to the 
cross section. In those models it is stringently required to include the bottom-quark contribution]^ 
The gluon fusion cross section is known at the next-to-leading order (NLO) in the SM including 
top- and bottom-mass effects 18,19 , in the 2HDM and in the MSSM including contributions from 
squarks and gluinos 20 -34 . For the top quark, an effective field theory approach can be applied 


in which the top quark is considered to be infinitely heavy and can be integrated out from the 
full theory. In this approximation, Higgs production has been calculated up to the next-to-NLO 
(NNLO) inclusively (^-37 as well as fully-differential (38 -40 


Electro-weak contributions and 
effects beyond NNLO in the heavy-top approximation have been studied in Refs. |41|-(50| for 


example, while there was a large effort 51,^ to push the accuracy to next-to-NNLO (N^LO) 


which has been succeeded very recently 53 


for both the inclusive cross section at the NNLO 54-59 and differential quantities 60,61 


Finite top-mass effects have been shown to be small 

as long 


as no kinematical scale (such as the transverse momentum of a particle) that is not integrated 
out exceeds the top-mass threshold. 


The full dependence of the top- and the bottom-mass at the NLO has been included so far 
in a POWHEG-type |62| matching to parton showers (PSs) [33] , the analytically resummed 
transverse momentum spectrum of the Higgs boson at NLO-I-NLL 63 , a MC@NLO-type (641 


'^See Refs. for a theoretical overview. 

^Note that m theories with an enhanced bottom Yukawa the associated Higgs production with bottom quarks 
becomes relevant, see Refs. [9 17 and references therein. 
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matching to the Herwig Monte Carlos [65 -67 
differential NNLO 


approach [70,71 


the NNLO+NNLL jet-vetoed 68 and the fully- 


[69] cross section; and in some approximated form recently also in the NNLOPS 
Furthermore, the 2HDM as well as supersymmetric effects from squarks and 
gluinos within the MSSM [20f|34[ have been implemented in the first two approaches from that 
list [33y72[. In this manuscript, we report on a new implementation of NLO QCD corrections in the 


SM, 2HDM and MSSM applying the MC@NLO-type matching to both Herwig and Pythia showers. 
We work in the framework of MadGraph5_aMC@NL0 73 and combine its capabilities with the 
corresponding amplitudes provided by SusHi [74] . The linking of SusHi to MadGraph5_aMC@NL0 


is realized by a scripl|^ aMCSusHi. Its usage as well as the application of the combined code to 
obtain cross section predictions in the SM, the 2HDM and the MSSM is detailed in this paper|^ 


The manuscript is organized as follows: In Section]^ we present a brief overview of the elements 
of the computation at hand. Section is dedicated to introduce aMCSusHi and is separated 
in three parts which cover: how to use the script (Section |3.1[ ), how to run the resulting code 
(Section |3.2| ) and how to treat the shower scale (Section |3.3| ). We will show a brief application of 
the code to phenomenological results in Section]^ and conclude in Section]^ 


2 Outline of the calculation 


The goal of this paper is to present a tool which allows for the computation of arbitrary infra-red 
safe differential observables at both the parton- and the hadron-level for the production of neutral 
Higgs bosons via gluon fusion in the SM, the 2HDM and the MSSM by matching the NLO cross 
section to a shower. 


The relevant NLO matrix elements are taken from Ref. 74 , which include both SM-like contri 


butions and sbottom, stop and gluino effects. Examples of corresponding Feynman diagrams 
are illustrated in Fig.[^ They are combined and matched to a parton shower by the well- 
known MC@NLO-method. The matched cross section in the MC@NLO framework can be written 
symbolically as: 


®‘-^/MC@NLO J 


Bn + Vn + j 


zno) 


+ / [d^n+lRn+1 - -^nUO) , 


( 1 ) 


where Bn determines the Born-level cross section, W the virtual (including mass factorization) 
and Rn+i the real corrections; is the Monte Carlo subtraction term, with the same 

IR poles as Rn+i, the Monte Carlo phase space d^^^^ tends to d4)„_|_i in the IR limits, and 
^(^Mc _ The quantity is the shower spectrum for observable O, as obtained 

by running the shower starting from an n-body configuration. 

^aMCSusHi can be downloaded under https: //cp3. irmp.ucl.ac .be/pro -jec ts/madRraph/wiki/HSushi#nol 
■^SusHi has recently been extended to include NMSSM Higgs production [^, which may be made available in 
aMCSusHi in the future. 
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(g) (h) (i) 

Figure 1: A sample of Feynman diagrams for gg ^ (j) contributing to the NLO cross 
section; (a-c) LO, (d-g) virtual and (h-i) real corrections. The graphical notation for 
the lines is: solid straight A quark; spiraled = gluon; dashed = scalar (squark or Higgs); 
spiraled with line = gluino. 


The cross section at Born-level is derived from the LO diagrams for gg ^ (j) where 4> € {h, H, A}, 
see e.g. Fig.[^(a)-(c). The NLO virtual and real corrections are governed by diagrams like the 
ones shown in Fig.[^(d)-(g) and Fig.[^(h)-(i), respectively, and similar ones with quark loops 
replaced by squark loops. 


Eq. Q is implemented for all standard parton showers 76-81 in the fully-automated framework 
MadGraph5_aMC@NL0. This code determines NLO QCD corrections to arbitrary scattering processes 
at the LHC. On the basis of UFO models [82| , the code even allows to carry out computations 
in any theory beyond the SM in a general manner as soon as the renormalization is known and 
implemented in a UFO model, see Refs. 73,83-85 for further information. However, the Higgs 
production mode through gluon fusion is a special one being loop-induced already at the LO. Such 
processes cannot be handled in a fully-automated manner by any code to date, since it requires 
the automation of two-loop amplitudes which is beyond current technology. Therefore, we have 
treated Higgs production through gluon fusion in the SM, 2HDM and MSSM as a special case, 
by linking the relevant amplitudes from SusHi. Furthermore, as far as the MSSM is concerned 
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SusHi requires a link to FeynHiggs 86 -96 which evaluates the corresponding Higgs masses and 
couplings in user dehned scenarios. Setting up the SusHi amplitudes in MadGraph5_aMC@NLD is 
handled by a publicly available script called aMCSusHi, which is automated to create the gg ^ (j) 
process folder; download SusHi and FeynHiggs; compile, install and link them; and replace the 
relevant amplitudes in the process folder. In the upcoming section, we describe the application 
of the script and explain the necessary steps to obtain phenomenological results. 


3 aMCSusHi script 

This code is based on MadGraph5_aMC@NL0 and SusHi. For further information on theses codes 
we refer the reader to corresponding publications (73|[74| . After using the script to set up the 
code, we will focus on the relevant user inputs to obtain phenomenological predictions for Higgs 
cross sections in the SM, the 2HDM and the MSSM. 

3.1 Usage of the script 

The aMCSusHi script is available for download from the website https : //cp3. irmp. ucl. ac. be/ 
projects/madgraph/wiki/HSushi#nol, It is fully automatic in setting up the gg ^ 4> process 
folder which includes downloading, installing and linking FeynHiggs and SusHi. At first, a 
dummy HEFT process folder for gg ^ h m. the five-flavor scheme at NLO (without the virtuals) 
is created^ then, the HEFT amplitudes in MadGraph5_aMC@NL0 are consistently replaced by the 
ones from SusHi (including the virtuals). The set-up requires only a single call of the aMCSusHi 
script: 

> ./set_up_ggH_MSSM_script.pl <ggH-folder> [<FeynHiggs-folder> [<SusHi-folder>]] 


The hrst argument is mandatory and determines the path to the process folder for gg ^ cf) that is 
generated by the script. The only requirement is that this folder has to be dehned as a sub-folder 
of the MadGraph5_aMC@NL0 directory. The second and third arguments are optional. If there 
are compiled versions of FeynHiggs and SusHi available on your computer you can give the 
names of the folders that contain the hies libFH.a and libsushi.a, respectively. When executed 
with two (one) arguments the script will ask whether it should automatically download and 
install SusHi (and FeynHiggs). The script will always download the latest versions of these 
codes. While running, the script requires some user inputs: It asks whether or not SusHi (and 
FeynHiggs) should be downloaded, in which folder they should be installed (default is inside the 
<ggH-folder>) or where to hnd SusHi (and FeynHiggs) if already installed. The user is simply 
required to follow these on-screen instructions. Furthermore, the script creates log-hles in the 
working directory for the download (” XX_curl.log “), the configure command (” XX_conf.log “) 
and the compilation (”XX_make.log“), where XX=FH for FeynHiggs and XX=SusHi for SusHi. 

®See the README file how this can be easily obtained with the MadGraph5_aMC@NL0 script. 
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These files are supposed to give complementary information for any kind of troubleshooting]^ 
Further information about the aMCSusHi script provides the README file. 

3.2 Running the code 

Once the gg —)■ 4> process folder has been set up by the script, the run can be started directly 
from the <ggH-folder> by typing 

*■ "V 

> ./bin/generate_events 

and following the usual steps in MadGraph5_aMC@NL0 to choose the run-modes (order, shower or 
fixed-order, madspin). Before running the code one may want to modify the input settings. In the 
following we will discuss the differences between aMCSusHi and the ordinary MadGraph5_aMC@NL0 
code regarding the input files. 

The input cards can be found under <ggH-folder>/Cards/, where the param_card.dat, run_card.dat 
and shower_card.dat contain all the essential information. The run card controls the usual pa¬ 
rameters, e.g. the renormalization (/xr) and factorization scale (/xp). Note, that by default 
the flags f ixed_ren_scale and fixed.fac.scale are set to false, so that these two scales 
are chosen on an event-wise basis. Their values are specified in the last routine of <ggH- 
folder>/SubProcesses/setscales.f which, due to the default option dynamical.scale.choice = —1 
in the run card, sets /xr = /xp = Ht/2 = 1/2 where i runs over all final state 

particles and m* and pr(x) are their mass and transverse momentum, respectively. This choice is 
reasonable, since it respects effects from hard radiation and corresponds to a value of in 

the soft/collinear limit which is the current recommendation for the total inclusive gg ^ (j) cross 
section [^[^. 

Also the shower card in aMCSusHi contains no new information and has the usual functionality. 
Since MadGraph5_aMC@ML0 supports all standard parton showers, for the first time PythiaG and 
PythiaS can be applied at NLO+PS to SM Higgs production in the full theory in the MC@NLO 
framework. In general, it is advisable to apply the most recent versions of the showers for 
meaningful physics runs. 

The param_card.dat, on the other hand, receives some significant changes by the aMCSusHi 
script. The new version basically combines the orignial parameter card from MadGraph5_aMC@NL0 
with the input hie from SusHi which are both written in the SUSY Les Houches accord (SLHA) 
format |97| and, therefore, easily connectable. We will address the different options in the 
param_card.dat in more detail, since there are a number of changes and some of the original 
parameters loose their functionality. In the SLHA format inputs are organized in blocks which 
have different entries that are characterized by a number. For simplicity, we introduce the 
following short-hand notation: Block example [i] corresponds to entry i in Block example. 

®Please note that even in cases where the compilation of FeynHiggs or SusHi fails, the linking may still work 
fine as long as the files libFH.a and libsushi.a have been created in the corresponding library folders. 
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E.g., entry 25 of Block mass (Block mass [25]) in the SLHA format is devoted to the Higgs 
mass in the SM, which is required as an input in the param_card.dat. A typical parameter card 
of aMCSusHi in the SM is shown below: 


################################### 

## INFORMATION FOR MASS 
################################### 

Block mass 

15 1.7770006+00 # MTA 

23 9.1188006+01 # MZ 

25 1.2500006+02 # MH -- only eff6Ctiv6 if FEYNHIGGS Block is absent 

[. . .] 

################################### 

## INFORMATION FOR SMINPUTS 
################################### 

Block sminputs 

1 1.3250706+02 # aEWMl 

2 1.166390e-05 # Gf 

3 1.1800006-01 # aS 

# additional information needed for SusHi 

4 9.1187606+01 # m_Z(pole) 

5 0.416000e+01 # m_b(m_b) -- only used if m_b is not on-shell 

6 1.730000e+02 # m_t(pole) -- top mass is set here 

################################### 

## INFORMATION FOR YUKAWA 
################################### 

Block yukawa 

15 1.777000e+00 # ymtau 

################################### 

## INFORMATION FOR DECAY 
################################### 

DECAY 6 1.491500e+00 # WT 

DECAY 23 2.441404e+00 # WZ 

DECAY 24 2.047600e+00 # WW 

DECAY 25 6.382339e-03 # WH 

DECAY 9000006 6.382339e-03 # WHl 

[. . .] 


Higgs (H) 

2 = MSbar(muR) 


################################### 
## INFORMATION FOR SUSHI 
################################### 
Block sushi 


1 

2 

Block 

1 

4 

Block 

1 

2 

3 


0 

0 

renormbot 

0 

4.75d0 
factors 
0. dO 
1. dO 
1. dO 


model: 0 = SM, 1 = MSSM, 2 = 2HDM 

0 = light Higgs (h), 1 = pseudoscalar (A), 2 = heavy 

Renormalization of the bottom sector 


m_b used for 
mbOS fixed -- 


bottom Yukawa: 
used if m_b is 


0 = OS, 1 = MSbar(m_b) 
on-shell (default) 


factor for yukawa-couplings: 

t 

b 


Most of the inputs are self-explanatory due to the comments initiated by the hash symbol # after 
the entries. Furthermore, the standard SLHA blocks match the universal convention of Ref. [97]. 
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Some of the inputs, though, require further comments. In the Block mass all parameters have 
the expected function, except for the top and the bottom mass. Block mass [6] and Block 
mass [5], respectively. While the former only affects and is required for the shower, the latter 
can be omitted completely. Instead, due to the link to SusHi the top mass that is used for the 
top loop and the top Yukawa is set in Block sminputs [6] and is expected to be on-shell. For 
the bottom mass, on the other hand, SusHi allows for three different choices: on-shell scheme or 
MS scheme with mfe(mft) or mb(/iR), which can be switched in Block renormbot [1] by a value 
between 1 and 30 Also here the recommendation is to use the on-shell scheme which according 
to Refs. [^,98 assures the cancellation of large logarithms ln(m/i/mb) at NLO QCD, while the 
MS scheme does not, due to an incomplete resummation of these terms. The on-shell mh value is 
determined by Block renormbot [4], while when the MS scheme is chosen the corresponding 
input of mb{mb) is set in Block sminputs [5]. The other entries of Block sminputs again have 
the same impact as in the usual MadGraph5_aMC@ML0 code. The same is true for Block yukawa. 
For the decay of light and heavy Higgs bosons one may specify a finite width of the Higgs boson 
in the respective BSM scenario by using Decay 25 irrespective of wether the light or the heavy 
Higgs boson is considered. At this point we shall remark that the particle identification number 
(PID) in the generated event files is always 25 regardless of the Higgs boson under consideration]^ 
This is irrelevant for the production (which is correctly computed through the SusHi amplitudes), 
but plays a role for the decay where the shower will consider particle 25 to be the light Higgs, 
which is indeed fine for any scalar Higgs, but a problem for pseudo-scalar ones. Therefore, decays 
of a pseudo-scalar Higgs are currently not supported in the official version of aMCSusHi. A user 
interested in decaying the pseudo-scalar Higgs is strongly encouraged to contact us. 


The other parameters are relevant to SusHi. Block sushi [1] chooses the model with the three 
options SM (Block sushi [1] = 0), MSSM (Block sushi [1] = 1) and 2HDM (Block sushi [1] = 
2). The second entry of Block sushi determines the Higgs boson: light Higgs (Block sushi [2] = 
11 or 0), pseudo-scalar Higgs (Block sushi [2] = 21 or 1) and heavy Higgs (Block sushi [2] = 12 
or 2). The choice of the masses of the relevant Higgs bosons depends on the model. As stated 
before, Block mass [25] sets the Higgs mass in the SM. In the 2HDM, this entry corresponds 
to the mass of the light Higgs boson, while Block mass [35] and Block mass [36] specify the 
input for the heavy and the pseudo-scalar Higgs, respectively. All other 2HDM inputs are set in 
the information for SusHi. For reference, we give an example of the corresponding inputs for a 
heavy Higgs in the 2HDM below0 

[. . .] 

################################### 

## INFORMATION FOR SUSHI 
################################### 

Block sushi 

12 # model: 0 = SM, 1 = MSSM, 

22 # 0 = light Higgs (h), 1 = 


2 = 2HDM 

pseudoscalar (A), 2 = heavy Higgs (H) 


^Note that only in the on-shell scheme the automatic MadGraph5_aMC@NL0 reweighting for is functionable. 
®Bear in mind that this has to taken into account to identify the Higgs in the analysis of the showered events. 


®A link to 2HDMC 
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with the corresponding input convention is currently not supported. 










Block renormbot # Renormalization of the bottom sector 

10 # m_b used for bottom Yukawa: 0 = OS , 1 = 

22 # tan(beta)-res. of Y_b: 0 = no, 1 = naive 

4 4.75d0 # mbOS fixed 

Block 2hdm # 2HDM version according to arXiv:1106.0034 

2 # (1=1,2=11 ,3=III ,4=IV) 

Block minpar 

3 50d0 # tanb 

Block alpha 

0.0247d0 # mixing in Higgs sector 

Block factors 


1 OdO 

2 IdO 

3 IdO 


# factor for yukawa-couplings : c 

# t 

# b 


MSbar(m_b), 2 
2 = full (for 


MSbar(muR) 
OS only) 


Additionally to the inputs which we defined already for the SM the following parameters have to 
be set in the 2HDM: Block renormbot [2] specifies whether or not a resummatior of terms 


enhanced by tan(/3) is applied through reweighting of the bottom Yukawa coupling as described 
in Ref. |74| ; Block 2hdm determines which type of the 2HDM is used; the value of tan(/?) is set 
through Block minpar [3]; and the mixing angle a corresponds to the entry in Block alpha. 


The computation of MSSM Higgs cross sections requires Block extpar, Block f eynhiggs and 
Block renormsbot in addition, which fix the parameters of the third family of quarks and 
squarks, determine the FeynHiggs inputs and yield information on the renormalization of the 
sbottom section, respectively. We will not provide any further information on these blocks, 
instead, we refer to the SusHi manual (Tog) and the FeynHiggs man pages dzl- Moreover, 
Block alpha can be omitted in the MSSM and the Higgs masses in Block mass have no effect, 
since they are determined by FeynHiggs, once Block feynhiggs is present. The MSSM Higgs 
mass that has been computed and applied in the run is provided to the user in Block mass [25] 
of the parameter card, which will be overwritten by the mass of the respective Higgs boson at 
the beginning of each MSSM run. 


So far we did not comment on the Block factors. It allows to turn on and off individual 
contributions in all models. In fact, it even provides the possibility to rescale the respective 
Yukawa couplings by choosing values different from zero and one. With Block factors [1] 
one can include the charm quark in the computation. This requires to specify its MS mass 
mc{mc) in Block sminputs [8] which is then translated to its on-shell mass. Furthermore, Block 
factors [2] and Block factors [3] multiply the top and the bottom Yukawa, respectively. In 
the MSSM, the stop Yukawa is rescaled by Block factors [4] and the sbottom one by Block 
factors[5]. 


For further information on the input cards we refer the reader to Ref. 73 of MadGraph5_aMC@NL0 


the manual of SusHi [Tog] and the man pages of FeynHiggs |107| . Three example parameter 
cards are provided in the folder <ggH-folder>/Cards; one for the SM (param_card.dat_SM), the 
2HDM (param_card.dat_2HDM_scenB) and the MSSM (param_card.dat_MSSM_mhmodp). They 
match the scenarios that we study in the result section of this paper. 


i°See Refs. 
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for further information. 
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3.3 Choosing different shower scales 


The choice of the shower scale is a very peculiar one in the gluon fusion process. In presence 
of the bottom-quark loop, factorization of soft and collinear radiation maybe spoiled at scales 
significantly smaller than the Higgs boson mass. This was pointed out by Ref. in the context 
of analytic transverse momentum resummation. On the other hand, these terms might well be 
treated as a hnite remainder as long as their impact remains moderate [68| . 

Due to their additive matching of the resnmmed low-p^ region with the fixed-order distribution 
valid at large transverse momenta, analytic p^-resummation and the MC@NLO method are quite 
similar. In both cases there is a scale associated with that matching, the resnmmation scale 
Qres and the shower scale Qshi respectively. These scales can be interpreted as transition scales 
that separate the soft/collinear from the hard physics, very similar to the factorization scale of 
the PDFs. In other words, they define the range where resummation, and therefore the shower 
in MC@NLO, takes effect. Their value has to be chosen of the order of the typical scale of the 
problem. 


In gluon fusion, the typical scale depends on the quark considered in the loop. Since mt ~ m^, 
there exist only two relevant scales for the top-quark loop {rUfj) and pt) and the shower scale can 
be chosen of the order of the Higgs mass. When considering the bottom loop, on the other hand, 
we face a three-scale problem {m^, mi, and pt) which has not been solved to date. However, it 
has been suggested 69 to apply a lower transition scale to the bottom contribution, which, in 


particular, respects the fact that soft/collinear factorization is valid only up to smaller scales for 
the bottom loop. In Ref. 72 it was further proposed to separate three contributions according 


to their Yukawa couplings: the square of the top and the bottom, and their interference; and 
choose separate shower/resummation scales for all of them. This splitting allows for a model 
independent treatment of the problem by a rescaling of the individual contributions with the 
respective top and bottom Yukawas of a specihc scenario in the 2HDM as well as the MSSM when 
neglecting squark effects. In the literature, two pragmatic approaches with physical motivation 
have been presented [^,108 to determine separate scales for the three contributions. Their 
comparison will be studied elsewhe re [109] . When studying phenomenological results in Section 
we will apply the scales from Ref. [72[ (referred to as ”HMW“ in what follows). 


The separation of the bottom contribution (including the interference) from the top one with 
different shower scales {Qh and Qt, respectively) requires three runs in MadGraph5_aMC@NL0, 
which have to be combined as follows: 


Qh) = MQt) + (^t+b{Qb) - MQb)- (2) 

To obtain all three contributions of different Yukawa origin with different scales that allows for a 
model independent treatment, on the other hand, MadGraph5_aMC@NL0 has to be run hve times: 

Qbj Qth) — T ^biQb') T (^t+b{Qtb') ^t{Qtb) ^b{Qtb) ■ (3) 

The scales Qt, Qb and Qtb determine the scale for top, the bottom and their interference. 
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respectively. As indicated before, the individual contributions can be separated using the Block 
factors in the parameter card. 

The shower scale in MadGraph5_aMC@NL0 cannot simply be accessed through the input cards, 
since it requires an advanced user to be familiar with its specific treatment in the code. 
MadGraph5_aMC@NLD does not use a simple fixed scale for Qshj instead, it statistically extracts the 
shower scale from a distribution peaked at a specific value. The user can only change the range 
of the interval of the distribution which of course also affects the peak. Therefore, we identify 
Qt, Qb and Qth in Eqs. © and (§ with the peak of the respective shower scale distributions. 

The so-called shape parameters define the interval of the distribution from which the shower scale is 
picked on a event-wise basis. They can be specified in the include file <ggH-folder>/SubProcesses/ 
madfks_mcatnlo.inc, where the relevant part is given by (MadGraph5_aMC@NLD default values): 

[. . .] 

c Define lower and upper veto range (see MC subtraction terms) 
double precision frac_low,frac_upp 
parameter (frac_low=0.IdO) 
parameter (frac_upp=l.OdO) 

c Other control switches on veto (see MC subtraction terms) 
double precision scaleMClow,scaleMCdelta 
parameter (scaleMClow=10.dO) 
parameter (scaleMCdelta=20.dO) 

[. . .] 


The parameters frac_low and fracmpp are used to compute the upper and lower bounds of the 
Qsh distribution, which will be explained in more detail below, while scaleMClow allows to set 
an absolute value of the lower bound on Qsh and scaleMCdelta is used to apply a minimal value 
to the size of the distribution interval. In formulas the interval is defined bfH 

Qrain ^ Qsh ^ Qmaxj with 

Qmin = max(frac_low • ^/sQ, scaleMClow) and (4) 

Qmax = max(frac_upp • Qmin + scaleMCdelta), 


where sq is the Born-level partonic center of mass energy squared. Evidently, scaleMClow and 
scaleMCdelta only take effect if the interval obtained through frac_low and frac.upp does not 
meet the corresponding restrictions. The corresponding Qsh distribution is peaked around 


Qpeak ~ 


(frac_low + fracmpp) —y 


(5) 


For a 2 —)■ 1 process like gluon fusion this relation is an identity and equals the mass of the 
final state particle, i. e., the Higgs mass in the case of gluon fusion. To change the peak-value 
to its half, e. g., for shower scale variations, one can simply divide frac.low and f rac.upp by a 


^'^For further information we refer the reader to Section 2.4.4 of Ref. 


73 
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factor of two. In this sense, it is convenient to keep the ratio between frac.low and frac.upp 
a constant, which in the default setup of MadGraph5_aMC@NLD is a factor of ten. Under this 
prerequisite, in order to choose a specific shower scale Qsh for gg ^ 4>we simply have to determine 


frac_upp 


2 Qsh 

1.1 


and frac.low 


frac.upp 
10 


( 6 ) 


Here and in what follows, we associate Qpeak with the shower scale Qsh and vice versa unless 
stated otherwise. After modifying the corresponding include file accordingly, MadGraph5_aMC@NL0 
has to be recompiled. This can be achieved by typing 


> make clean 

- 


inside the <ggH-folder> which forces a recompilation during the next run of the code. 

In Section we show some applications of the aMCSusHi code and study the effect of different 
treatments of the shower scales. 


4 Results: brief application 


The gg ^ (j) process folder created by the aMCSusHi script preserves all the highly convenient 
features that come with MadGraph5_aMC@NL0. Besides many others, this entails an interface to the 
most common showers, the fully-automatic determination of scale and PDF variations without 
any extra-costs of computing time 110 , the creation of any number of observables with a single 


run and analysis routines available for the most important processes including the gluon fusion 
Higgs production mode. 

aMCSusHi allows to compute gluon-induced Higgs production including the complete dependence 
on the quark masses in the SM for the first time in a MC@NLO-type matching applying all versions 
of Pythia and Herwig Monte Carlos. While previous computations did only feature the Herwig 
65 -67 , phenomenological results exist to our knowledge only for Herwig6 [111 


showers 


As 


first application we therefore study the impact of different showers on the top- and bottom-mass 
effects with respect to the heavy-top approximation at the 13 TeV LHC. For this purpose, Fig.[^ 
shows the ratio of the NLO-PPS computation including mass effects and the corresponding cross 
section in the heavy-top limit as a function of the transverse momentum of the Higgs for different 
Monte Carlos: PythiaS (black, solid), Herwig++ (red, dotted), Pythia6 pT-ordered (blue, dashed 
with points), PythiaS Q-ordered (green, dash-dotted with open boxes) and HerwigS (yellow, solid 
with filled boxes). We apply the MSTW2008 68%CL NLO PDF set |112] with the corresponding 
value of the strong coupling constant. The shower scale has been chosen as Qsh = in 

all cases, while for /rp and /ip we use the defaults specified in Section [T^ Clearly, the mass 
effects are hardly dependent on the specific Monte Carlo which is particularly evident at small 
{pt ^ 50 GeV) and large [px ^ 150 GeV) transverse momenta. Nevertheless, there are some 
visible differences in the intermediate region which consistently discriminate the Herwig from 
the Pythia showers. Overall, they are at most 5% though and therefore still moderate. 
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da/dpT(c|)) / dOht|/dpT((l)) 



Figure 2: Transverse momentum distribution of a SM Higgs at NLO+PS in the full 
theory normalized to the one in the heavy top effective field theory for different Monte 
Carlos: PythiaS (black solid), Herwig++ (red, dotted), Pythia6 p^-ordered (blue, 
dashed with points), Pythia6 Q-ordered (green, dash-dotted with open boxes) and 
Herwig6 (yellow, solid with filled boxes). 


Fig-i shows the effects of quark masses with respect to the heavy-top approximation as well, 
but for different choices of the associated shower scales. In all cases, the denominator and 
therefore the distribution in the heavy-top limit is computed with the respective scale of the 
top contribution Qsh = Qt- As we observed before the Monte Carlo dependence is quite small; 
therefore, we only consider the PythiaS shower. For reference the black solid curve is the same 
as in Fig.^ with Qsh = for all contributions. We compare it to the scales choices proposed 

in Ref. , which imply setting the shower scale of the bottom contribution (including the 
interference) to the bottom mass following Eq. ([^ (red dotted curve). For the blue dashed curve 
with points we chose the HMW scales determined in Ref. 72 which can be found in Table 1 of 


that paper, applying a three-scale approach according to Eq. ([^ (Q^ = 49 GeV, Qth = 34 GeV 
and Q\) = 23GeV). The green dash-dotted curve with open boxes serves mostly for comparison 
with previous Herwig6 results [m] which were computed with the scales of Ref. 


69 as well. 


For the pT distribution in Fig.[^(a), the change of the scale of the bottom contribution to Qh = rnt, 
has a significant impact on the mass effects at small and intermediate transverse momenta. It 
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Figure 3: (a) Same as Fig.[^ but for different choices of the shower scales, see text for 
details; (b) corresponding plot for the rapidity distribution of the Higgs. 


develops an extremely steep drop at small transverse momenta which due to unitarity affects 
also the intermediate pT-range in the opposite direction. The beneht of the usage of such a low 
scale is clearly disputable. While the Herwig6 curve agrees rather well with previous result of 
Ref. m becoming flat for pt ^ 5 GeV, the PythiaS curve develops a steep increase in this 
region. This signals a significant Monte Carlo dependence at very small pT which is not observed 
for larger Qi, scales. Furthermore, the rigorously low value also poses a technical problem in the 
code regarding the fact that the default shower scale choice in MadGraph5_aMC@NL0, as explained 
in Section is a distribution. In order to solve this problem, we had to use a fixed value of 
Qsh = Qb by setting frac_low = fracmpp = Qblrrih and scaleMClow = scaleMCdelta = 0. 

Considering the HMW scales, the scale of the bottom contribution is not chosen at such low 
values. We hnd that the mass effects in this case (blue dashed line with points) are rather similar 
to the ones where all scales are set to m/i/2 (black solid line), although the individual HMW 
scales being quite different from this value. Looking at the rapidity distribution in Fig.|^(b), on 
the other hand, we observe the expected feature of being essentially insensitive to any choice of 
the respective shower scales. We shall note at this point that simply due to their inclusion in the 
default analysis we were able to produce a large number of further observables at no additional 
computing cost. 

To demonstrate the range of applicability of aMCSusHi, we consider two realistic BSM scenarios in 
Fig. 0 and Fig.[^ the heavy Higgs boson in Scenario B of Ref. (HD (a bottom dominated 2HDM 
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Figure 4: Transverse momentum distribution of (a) the heavy Higgs boson and (b) the 
associated hardest jet computed in a bottom dominated scenario of the 2HDM (see text 
for details). Graphical notation is the following: black solid curve shows PythiaS at 
NLO+PS, red dotted curve is the same at LO+PS (normalized the NLO) and the blue 
dashed one with points corresponds to the fixed order curve at NLO. 


scenario); and the pseudo-scalar Higgs boson in the "^(800,40) MSSM scenario 
Table 2 of Ref. Wa- The corresponding input files can be found in the folder <ggH- 


114 


defined in 


blder>/Cards. 


In Fig.we study the transverse momentum distributions of the Higgs and the hardest jet, while 
Fig. [^depicts their rapidity distributions. In both cases we apply the HMW scales of Ref. (7^ . 
At low transverse momenta the pT distributions have similar shapes comparing the red dotted 
(LO-I-PS) to the black solid curves (NLO-l-PS). This can be easily inferred from the first inset 
where all curves are normalized to the black solid line in the main frame. However, it is well 
known that at the LO-I-PS pT distributions yield unphysical results for transverse momenta 
beyond the shower scales indicated by a steep drop. Note that both curves are normalized to 
the same (the NLO) cross section. Continuing the comparison at hand, we observe a significant 
reduction of the scale uncertainties shown in the lower inset, where the bands are obtained by 
dividing the upper and lower bound of the respective cross section by the same central cross 
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Figure 5: Rapidity distribution of (a) the pseudo-sc alar Higgs boson and (b) the 


associated hardest jet computed in the scenario 114 of the MSSM with Ma = 


800GeV und tan/3 = 40. Graphical notation is the same as in Fig.[^ All curves are 
normalized so that their bins add up to one. 


section as in the first inset. The uncertainties correspond to the independent variation of all 
unphysical scales (/tr, /ip, Qt-, Qb, Qtb) by a factor of two. Gomparing NLO-I-PS to the NLO 
fixed-order result denoted by fNLO, we observe the expected matching towards large transverse 
momenta. 


In order to compare shapes, the rapidity distributions in Fig.are normalized in a way that the 
sum of their bins yields one. We see that for the Higgs rapidity in Fig.[^(a) all curves agree 
extremely well in terms of shape up to the forward region in which, nevertheless, the deviations 
are still well within the respective uncertainty bands For the rapidity distribution of the 


hardest jet the same is true for the two showered results, while the fNLO distribution, on the 
other hand, agrees only in the central region |?/(ji)| < 3, but features a significant enhancement 
when the hardest jet is more forward. In this region the cross section will receive large effects of 
collinear radiation which renders the shower to yield the more reliable description. 


^Note that the Higgs reaches its kinematical limit slightly before \y{4>)\ ~ 3 already. 
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5 Conclusions 


In this article we presented the new tool aMCSusHi which is a link between MadGraph5_aMC@NL0 
and SusHi for the computation of Higgs cross sections in gluon fusion at hadron colliders. The 
code gives NLO+PS accurate results in the SM, 2HDM and MSSM. The inputs in the MSSM 
are conveniently controlled through a link to FeynHiggs. We discussed the specific treatment 
of the shower scale in MadGraph5_aMC@NL0 and pointed out its special role in the context of 
gluon-induced Higgs production. In the phenomenological part we study the impact of different 
shower scale choices on the mass effects in the SM. Furthermore, we studied results for 2HDM and 
MSSM Higgs production as an application of aMCSusHi. The aMCSusHi script can be downloaded 
from https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/HSushi#nol, 

As an outlook, one may improve the SM prediction for the Higgs production mode through 
gluon fusion by merging the NLO-I-PS cross section for gg ^ h m. the full theory with higher 
multiplicities computed in the heavy-top effective field theory. This is certainly beyond the scope 
of the present paper and is left for a future publication. 
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